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Abstract 

Taking the solar and the atmospheric neutrino experiments into account we discuss the lepton 
flavor violating processes, such as r — » py or p, — > e'y, in the minimal supersymmetric standard 
model with right-handed neutrinos (MSSMRN) and the supersymmetric SU(5) GUT with right- 
handed neutrinos [SU(5)RN]. The predicted branching ratio of /x — ► e'y in the MSSMRN with the 
Mikheyev-Smirnov-Wolfenstein (MSW) large angle solution is so large that it goes beyond the current 
experimental bound if the second-generation right-handed Majorana mass M„ 2 is greater than ~ 
10 13 (~ 10 14 )GeV for tan/3 = 30(3). When we take the MSW small angle solution, the p, -> ej rate is 
at most about 1/100 of that of the MSW large angle solution. The 'just so' solution implies 10 _J of 
that of the MSW large angle solution. Also, in the SU(5)RN the large p — > e'y rate naturally follows 
from the MSW large angle solution, and the predicted rate is beyond the current experimental bound 
if the typical right-handed Majorana mass M N is larger than ~ 10 13 (~ 10 14 )GeV for tan/3 = 30(3), 
similarly to the MSSMRN. We show the multimass insertion formulas and their applications to 
r — » p,y and p — > ey. 



1 Introduction 



Introduction of supersymmetry (SUSY) to the Standard Model (SM) is a solution for the naturalness 
problem on the radiative correction to the Higgs boson mass. The minimal supersymmctric standard 
model (MSSM) is considered as one of the most promising models beyond the Standard Model. Nowadays 
the signal of supersymmetry is being searched for by many experimental ways. 

Lepton flavor conservation, lepton number conservation in each generation, is an exact symmetry in 
the SM, however, it may be violated in the MSSM jy]. The SUSY breaking mass terms for sleptons have 
to be introduced phenomenologically. Then, the mass eigenstates for sleptons may be different from those 
for leptons. This leads to the lepton flavor violating (LFV) rare processes, such as n — ► ej, t — > /J.J, and 
so on. In fact, the experimental bounds on them have given a constraint on the slepton mass matrices. 

The structure of the SUSY breaking mass matrices for sleptons depends on the mechanism to generate 
the SUSY breaking terms in the MSSM. One of the interesting mechanisms is the minimal supergravity 
(SUGRA) scenario. Similar to the slepton masses, arbitrary SUSY breaking masses for squarks are 
also strongly constrained from the FCNC processes, such as K° — T(° mixing. In the minimal SUGRA 
scenario, the SUSY breaking masses for squarks, sleptons, and the Higgs bosons are expected to be given 
universally at the tree level, and we can escape from these phenomenological constraints. 

However, the universality of the SUSY breaking masses for the scalar bosons is not stable for the 
radiative correction. Especially, if the physics below the gravitational scale M grav (~ 10 18 GeV) has the 
LFV interaction, the interaction induces radiatively the LFV SUSY breaking masses for sleptons. Then, 
the LFV rare processes are sensitive to physics beyond the MSSM ||. 

Recently, the Super-Kamiokande experiment has given us a convincing result j^] that the atmospheric 
neutrino anomaly Q] comes from the neutrino oscillation. From the zenith-angle dependence of v e and 

fluxes following neutrino mass square difference and mixing angle are expected, 

A<, x - (Kr 3 -ICr 2 )eV 2 , 

sin 2 26 VftVX > 0.8. (1) 

From the negative result for v^-v^ oscillation in the CHOOZ experiment || it is natural to consider 
vx = v T from above results, and the tau neutrino mass is given as 

m„ x ~ (3 x I(T 2 - I x KT^eV, (2) 

provided mass hierarchy m VT 3> m Vll . 
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The simplest model to generate the neutrino masses is the seesaw mechanism Q . The neutrino mass 
Eq. (Q) leads to the right-handed neutrino masses below ~ (10 14 — 10 15 ) GeV, even if the Yukawa coupling 
constant for the tau neutrino mass is of the order of one. This means that a LFV interaction exists below 
the gravitational scale. Then it is expected that the LFV large mixing for sleptons between the second- 
and the third-generations is generated radiatively in the minimal SUGRA scenario, and that the LFV 
rare processes may occur with rates accessible by future experiments R, f§. In fact, the branching ratio 
of r — > in the MSSM with the right-handed neutrinos can reach the present experimental bound |j) . 

The solar neutrino deficit fllPf may also come from neutrino oscillation between v^-v^. The Mikheyev- 
Smirnov-Wolfenstein (MSW) solution due to the matter effect in the sun is natural for its explanation, 
and the observation favors 



Am 2 ei , y 



~ (8 x 1CT 6 - 3 x l(T 4 )eV 2 , 
sin 2 26 VeVY > 0.5, (3) 



Am 2 el/y ~(4x 10~ 6 - 1 x 10~ 5 )eV 2 , 

sin 2 20^ ~ (10" 3 -10- 2 ). (4) 

If the solar neutrino anomaly comes from so-called 'just so' solution neutrino oscillation in vacuum, 
the following mass square difference and mixing angle are expected [fl3|| , 

Ato 2 ^ ~(6x 10~ n - 1 x 10- 10 )eV 2 , 

sin 2 20 VtVY > 0.5. (5) 

It is natural to consider vy = u^, combined with the atmospheric neutrino observation. If one of the 
large angle solutions for the solar neutrino anomaly is true, the large mixing 6 V Ve may imply the LFV 
large mixing for sleptons between the first- and the second-generations. 

In this article we investigate the LFV processes in the supersymmetric models with the right-handed 
neutrinos, assuming the minimal SUGRA scenario. We take the above results for the atmospheric and 
solar neutrinos. In section 2 after discussing the origin of the observed mixing angles we calculate the 
t — > fij and jj, — > ej branching ratios under the assumption of the MSSM with the right-handed neutrinos. 
It is argued that the fj, — ► e'y rate depends on the solar neutrino solutions, and that especially the MSW 
large angle solution naturally leads to a large [i — > ej rate. In section 3 we consider them in the SU(5) 
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SUSY GUT with the right-handed neutrinos. Here also it is shown that the large \i — > rate naturally 
follows from the MSW large angle solution. Section 4 is for our conclusion. In Appendix A we give 
our convention used in this article. In Appendices B and C we show the multimass insertion formulas 
and their applications to r — > ^7 and /j, — > e-f, which are useful for estimating the LFV amplitudes 
and understanding the qualitative behavior of the LFV rates. In Appendix D the renormalization group 
equations (RGE's) relevant for our discussion are given. 

2 The Lepton Flavor Violation in the MSSM with Right-handed 
Neutrinos 

Before starting to investigate the LFV rare processes in the MSSM with the right-handed neutrinos 
(MSSMRN), we discuss the origin of the large mixing of neutrino between Vt-v^ or v^-v e . The MSSMRN 
is the simplest supersymmetric model to explain the neutrino masses, and following discussion is valid to 
the extension. The superpotential of the Higgs and lepton sector is given as 

IUmssmrn = U ij H 2 NiL j + f eij H 1 E i L j + -M ViVj N i N j +^H 1 H 2 , (6) 

where L is a chiral supcrfield for the left-handed lepton, and N and E are for the right-handed neutrino 
and the charged lepton. H\ and H 2 are for the Higgs doublets in the MSSM. Here, i and j are generation 
indices. After redefinition of the fields, the Yukawa coupling constants and the Majorana masses can be 
taken as 

feij = fe, $ij > 

M ViVj = U* k M Vk U\ p (7) 

where Vd and U are unitary matrices. In this model the mass matrix for the left-handed neutrinos (m„) 
becomes 

K)ij = V^ ik {rn v ) H V D ij, (8) 

where 

{m v )ij = m ViD [M~ 1 ]..m l , jD 

= VM ik m Vk V M kj- (9) 
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Here, m Vi r> — /^usin/3/%/2 and Vm is a unitary matrix.^] We assume f V3 ^ f V2 ^ f Vl , similar to the 
quark sector, and m Vr 3> m v 3> m Ue . Also, we take the Yukawa coupling and the Majorana masses for 
the right-handed neutrinos real for simplicity. 

When we consider only the tau and the mu neutrino masses, we parameterize two unitary matrices 

as 

/ cos9 D sin9 D \ ( cos 9 m sin 0m \ , lrV , 

V D = \ _. « „~„ , V M =[ _ fl .. • (10) 



sin#D cosOd J ' V — sin 0m cos9m 

The observed large angle 9 V Vt is a sum of 9d an d 9m- However, in order to derive large 9m we need 
to fine-tune the independent Yukawa coupling constants and the mass parameters. The neutrino mass 
matrix (m») in the second and the third generations is written explicitly as 



i/ 2 D ^u 2 D rn u 3 



(m) - - 1 = 1 (||) 



If the following relations are imposed, 



m l 3 p ^ m l 2 p ^ rn^pm^D 



(12) 



the neutrino mass hierarchy m Ur 3> m v and 9m — t/4 can be derived. However, it is difficult to 
explain the relation among the independent coupling constants and masses without some mechanism or 
symmetry. Also, the hierarchy m V3 r> 3> m V2 D suppresses the mixing angle 9m as 

tan2^2(^V^y (13) 



jn V3 D j \M V2V2 

as far as the Majorana masses for the right-handed neutrinos do not have stringent hierarchical structure 
as Eq. (jl^). Therefore, in the following discussion we assume that the large mixing angle between v T and 
comes from 9p and that U is a unit matrix. Similarly, it is natural to consider that the large mixing 
angle between and v e in the MSW solution or the 'just so' solution for the solar neutrino anomaly 
comes from Vd- 

The existence of the large mixing angles in Vd may lead to radiative generation of sizable LFV 
masses for the sleptons in the minimal SUGRA scenario. Though the SUSY breaking masses for the 
left-handed slepton are flavor-independent at tree level, the Yukawa interaction for the neutrino masses 
induces radiatively the LFV off-diagonal components in the left-handed slepton mass matrix. 

1 (hi) = (v cos j3/V2,0) T and (hi) = (0, v sin /3/v / 2) T with v ~ 246GeV. 
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The SUSY breaking terms for the Higgs and lepton sector in the MSSMRN are in general given as 

- £ SUSY breaking = {m\)ijV L J,Lj + {ml)ije Ri e Rj + (ml)ijV% i VRj 

+ml 1 h\hi +fh 2 h2 h\h 2 

+{A i jh 2 v R j Lj + K j hie* R J L j + ^v* Ri v Rj + B h hth 2 + h.c), (14) 

where II, e R , and v R represent the left-handed slepton, and the right-handed charged slepton, and the 
right-handed sneutrino. Also, h\ and h 2 are the doublet Higgs bosons. In the minimal SUGRA scenario 
at the gravitational scale the SUSY breaking masses for sleptons, squarks, and the Higgs bosons are 
universal, and the SUSY breaking parameters associated with the supersymmetric Yukawa couplings or 
masses (A or B parameters) are proportional to the Yukawa coupling constants or masses. Then, the 
SUSY breaking parameters in Eq. ( |l4| ) are given as 

(m|)y = (mf)ij = (m?) y = S tj ml, 

mil = ™H2 = m o> 
A u = fu i} a , A l i = f e ^a , 

S« = M VtVj bo, B h =(ib . (15) 

In order to know the values of the SUSY breaking parameters at the low energy, we have to include 
the radiative corrections to them. We can evaluate them by the RGE's. We present them in Appendix [d| 
and here we discuss only the qualitative behavior of the solution using the logarithmic approximation. 
The SUSY breaking masses of squarks, sleptons, and the Higgs bosons at the low energy are enhanced by 
gauge interactions, and the corrections are flavor- independent and proportional to square of the gaugino 
masses. On the other hand, Yukawa interactions reduce the SUSY breaking masses. If the Yukawa 
coupling is LFV, the radiative correction to the SUSY breaking parameters is LFV. The LFV off-diagonal 
components for (m|), (m|), and A l J are given at the low energy as 

«h - -^(3rnl + al)V£ ki V Dkj fllog^ 



(ml)ij ~ 0, 



3 „ ,„ „ „o , M s 



A% i - »/h,1/h ; ./;. log- : 



8n 2 — UKl o M 



where % ^ j. In these equations, the off-diagonal components of (rn|) and A e are generated radiatively 
while those of (m|) are not. This is because the right-handed leptons have only one kind of the Yukawa 
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K)32 - -g^(3mg + og)^33VD32/^ log (17) 



interaction / e and we can always take a basis where f e is diagonal. The magnitudes of the off-diagonal 
components of (m|) and A e are sensitive to f Vi and Vd-0 

As shown above, VD32 is expected to be of the order of one from the atmospheric neutrino observation. 
This leads to the non- vanishing (m|)32 and A^ 2 , which result in a finite r — > /i7 decay rate via diagrams 
involving them. The dominant contributions are proportional to 

As will be shown, if /„ 3 is of the order of one, the branching ratio of r — » /^7 may reach the present 
experimental bound. 

Moreover if Vo3i is finite, (m^jsi and (m~)2i are also large. They are approximately 
(m|) ai - - g ^(3mg+ag)yS,,V 1M1 ^Iog^, 

(m|) 21 * -^(amg+agjVSMV^Iog^. (18) 

This fact implies a sizable fj, —> ej rate because the amplitudes proportional to {rn 2 ^)23{m 2 z)3i or (m|)2i 
arc dominant. When Vd21 is also of the order of one to explain the solar neutrino anomaly, an extra 
contribution to (m£)2i nas to be taken into account as 

(m|) 21 ~ + a 2 ) (v£ 32 V mi f V3 log ^ + V* m2 V mi f V2 log ^) . (19) 

The experimental upper bound on the branching ratio of /i — > e7 is so severe that the predicted branching 
ratio may reach it even if f V2 is O(10 _1 ). 

2.1 The Branching Ratio of r — ► /i7 

Let us discuss the branching ratios of the LFV rare processes in the MSSMRN. First, r — > //7- The 
amplitude of the ef — > e^7 (i > j) takes a form 



ee Q *( g )^(p)»a a/3 /(Af )P L + A^P^ifr - q), (20) 



where p and q are momenta of and photon, and the decay rate is given by 



r«. <D = ^<(l4 J ' ) l 2 + l4 J ' ) l 2 )- c-'ii 



If (7 is not a unit matrix, the off-diagonal components for ("J^) and become 



(m f )ij ~ -—j(3m +a )V£ ki V Dl jf Vk f vl UZ m Ui m log 



O 9 ^ U UK I ^"■J""K,"^L Kill <■"" <J ji/T 3 

^ - -^^f^DMVDijUJ^U^U^log^L. (16) 



Then they are insensitive to the detail of U since the dependence on M Vi is logarithmic. 
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Here, we neglect the mass of ej. The amplitude is not invariant for the SU(2)i and U(l)y symmetry 
and the chiral symmetry of leptons. Then the coefficients A^f' and A^' are proportional to the charged 
lepton masses. Since in the MSSMRN the mismatch between the left-handed slepton and the charged 
lepton mass eigenstates is induced, A^ is much larger than A^p since A^ is suppressed by m ej /m ei 
compared with . Also, when tan/3(= V2/V1) is large, the contribution to A^ proportional to 
f e .V2(— — \[2m ei tan/3) becomes dominant. In the MSSMRN, the dominant contribution to r — ► ^7 is 
from the diagram of Figs. (|lj)(a) and (b) and its expression is 



A { ^ ] ~ m T ^/*M 2 taa/3(mf)82 

47T 



xD 



D 



1 2 



(M|,M 2 );m 2 



(™L™I T ), (22) 



which comes from the SU(2)z, interaction. The functions f n 2{x) and f C 2{x) are defined in Appendix |^ 
and the operator D[f(x); x] (x\, X2) to a function /(x) is defined by 

D[f(x);x](x u x2)= (/(xi) - /(as 2 )). (23) 

Here, for a demonstrational purpose, we take a limit where the SUSY breaking scale is much larger 
than the W and Z gauge boson masses and tan/3^ 1. This equation can be derived from the mass- 
insertion formula represented in Appendix ^|. The LFV A term can not give a dominant contribution 
when tan/3^ 1. 

In Fig. (||) we show the branching ratio of r — > fjr/ as a function of the left-handed selectron mass 
(mg L ). Here, m UT = 0.07eV, Vd33 = Vd21 = — Vd32 = Vd23 = 1/V2, and we assume that /„ 3 is as 
large as the Yukawa coupling constant for the top quark at the gravitational scale. This corresponds to 
M„ 3 ~ 10 14 GeV. Also, we impose the radiative breaking condition of the SU(2)ixU(l)y gauge symmetry 
with tan/3 = 3, 10,30 and the Higgsino mass parameter fj, positive. In our calculation we considered the 
experimental constraints from the negative results of the SUSY particle search. Though we do not assume 
the GUT's, we take the wino mass (M2) 130GeV and determine the other gaugino masses by the GUT 
relation for the gaugino masses. We use the formula for r — > /17 in Ref. j^] for the numerical calculation. 

The branching ratio is reduced where the left-handed selectron mass is comparable to the wino mass. 
This is because the slepton masses are almost determined by the radiative correction from the gaugino 
masses, and mg, which (to|)32 is proportional to, is negligible in the region. As mentioned above, the 



branching ratio is proportional to tan 2 j3 (see Eq. (p2[)), and the line for tan/3 = 30 is close to the 
experimental bound, Br(r — > ^7) < 3.0 x 10~ 6 Q. 
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In Fig. (0) we present the dependence of the branching ratio of r — > /17 on M V3 . Here, we take 



170GeV, and the other SUSY breaking parameters are the same as in Fig. 



The branching 



ratio is proportional to M% 3 since we fix m UT — 0.07cV. If 10 8 can be reached in the future experiments, 
we can probe M Ua > 10 13 (10 14 )GeV for tan/3 = 30(3).R 



2.2 The Branching Ratio of jx 



e7 



Next, we discuss /1 — » ej in the MSSMRN. The forms of the amplitude and the event rate are the same 
as those of r — > fi'y (Eqs. (^0 21)). This process has two types of the contribution, depending on the 
structure of the Yukawa coupling for the neutrino masses. One is the diagrams where (ml) 21 or A 2 , 1 is 
inserted, and another is those that (m'~)32 or A 32 and (m|)i3 or A)? are inserted. Then the dominant 
contributions (Figs. (||)(a)-(d)) are following, 



-m u — M2U tan 8 
A-k 



xD 



(m~ ) 2 i + 



(mf) 2 3(mf) 3 i 



(m|) 2 3(m 2 -) 3 




(M|,/?) (24) 



(m|-m| T ) 2 

Here, we take a limit where the SUSY breaking scale is much larger than the W and Z gauge boson masses 
and tan/3^ 1, again. We also assumed the mass degeneracy between the first- and the second-generation 
left-handed sleptons as 

2 2 2 2 _ 2 /r ) r\ 

m e L = ™Al = m £>e = TO iv = m £" ( 25 ) 

The functions f C 2,n2(x) and g C 2,n2(x) are defined in Appendix ^|. 

As mentioned above, if the solar neutrino anomaly comes from the MSW effect or the vacuum oscilla- 
tion with the large angle, Vd21 is expected to be large. This leads to non-vanishing (m|)2i- In Fig. (||), 
under the condition that 



0.91 0.35 0.24 
V D = ( -0.42 0.72 0.55 
-0.60 0.80 



(26) 



and = 0.004eV |17j we show the branching ratio of /1 — > ej as a function of M„ 2 . This corresponds to 
the MSW solution with the large mixing. Here we take Vd31 = and we will discuss a case with finite Vd31 
later. The input parameters are taken to be the same as in Fig. (||). For tan 3 = 30(3), the branching ratio 
reaches the experimental bound (Br(/i -> e 7 ) < 4.9 x lO" 11 El) when M V2 ~ 8 x 10 12 (8 x 10 13 )GeV. This 



An alternative way to prove (mi)32 is to search for the slepton oscillation |1EJ |lq ] . 
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corresponds to /„ 2 ~ 0.03(0.11). Future experiments are expected to reach 10 14 fig] . This corresponds 
to M v% ~ 10 11 (10 12 )GeV. 

If the solar neutrino anomaly comes from the MSW solution with the small mixing, we cannot dis- 
tinguish whether the mixing comes from Vd or Vm- If it comes from Vd, the branching ratio is smaller 
by about 1/100 compared with that in the MSW solution with the large mixing, as shown in Fig. (^). In 
Fig. (H) we assume that 

/ 1 0.04 0.03 \ 
V D = -0.04 0.79 0.59 (27) 
\ -0.60 0.80 / 

and m v = 0.0022eV fl7|| . Other input parameters are the same as Fig. (||). 
In Fig. §) we take 



V D = 



1 V2 2 2 
1_ 1 1 

V2 2 2 

V 3- 



(28) 



and m v = 1.0 x 10~ 5 cV [[l9). Other parameters are the same as in Fig. (||). This corresponds to the 
'just so' solution for the solar neutrino anomaly. Since the mu neutrino mass is smaller, the branching 
ratio is suppressed by 10~ 5 compared with that in the MSW solution with the large mixing. 

Next we discuss the branching ratio of \i — > ej when Vd31 is finite. In Figs. (||,^]) we show the branching 
ratio as a function of Vd31 and M V3 for tan /3 = 3 and 30. Here we assume that is negligibly small. 
The other parameters are the same as in Fig. (|^). The branching ratio is almost proportional to V^ 31 M^ 3 . 
Compared with this figure to Fig. (0), when M U2 — M U3 , the contribution from Vdsi 1S negligible in the 
MSW solution with the large mixing angle unless Vd3i is larger than 10 -2 5 . On the other hand, it can 
be dominant in the 'just so' solution even if Vd3i ~ 10 -4 . 

Finally we consider the /i + — > e + e _ e + process and the /i-e conversion on ||Ti. For these processes 
the penguin type diagrams dominate over the others, so the behavior of the decay rate is similar to that 
of /i — > e"f. For the /i —> 3e process the following approximate relation holds between the branching ratios 
of the two processes, 

a 8 / ml 11 \ 
Br( M ->3e) ~ -| (log^ - - J Br( M e 7 ) (29) 

~ 7 x 10~ 3 Br(^ -> e 7 ). (30) 
For the /x-e conversion rate T(/x — * e) a similar relation holds at tan/3 > 1 region, 

I> -> e) ~ 16a 4 Z 4 ff Z|F( (Z 2 )| 2 Br( A1 -» e 7 ). (31) 
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Here Z is the proton number in the nucleus, and -Zeff is the effective charge, F(q 2 ) the nuclear form factor 
at the momentum transfer q. For fl^i, Z e g — 17.6 and F(q 2 ~ — m^) ~ 0.54 ^l). We express the 
magnitude of the fi-e conversion with the normalization the muon capture rate in Ti nucleus. Then the 
normalized conversion rate R(fi~ — > e^jl^Ti) is approximately 

R(fi~ -> e-;||Ti) ~6x 10" 3 Br(/i -> &y). (32) 

The future experiment for the fi-e conversion is planed to reach R([f~ —> e~; |fTi) < 10~ 18 |2^| . 

3 The Lepton Flavor Violation in the SU(5) SUSY GUT with 
Right-handed Neutrinos 

In the SUSY GUT the gauge coupling unification is predicted, and the predicted weak mixing angle is 
consistent with the experimental data at the 1 % level of accuracy. Moreover if the unified gauge group 
is SO(10), the right-handed neutrinos are introduced automatically into the matter multiplet. However, 
in order to accommodate the observed large mixing angle in the framework of the SO(10) SUSY GUT, 
one needs unnatural extension of the simplest version of the SO(10) SUSY GUT. Hence in this article 
we do not discuss the SO(10) SUSY GUT. Here we investigate the SU(5) SUSY GUT with the right- 
handed neutrinos as one of the extension of the MSSMRN in which the small neutrino mass is naturally 
obtained and the large neutrino mixing angle is possible without unnatural fine-tuning. We here call this 
model as SU(5)RN, for brevity. After introducing the model we estimate the off-diagonal elements of the 
slepton soft mass matrices using the one-loop level RGE's under an assumption of the minimal SUGRA 
scenario. With them we study the LFV processes r — ► prf and /i — > ej. After that we comment on the 
b — > S7 branching ratio. We show that the LFV rates in this model is larger in general than those in the 
MSSMRN model, due to the fact that in this model the right-handed slepton mass matrix also can have 
non- negligible off-diagonal elements, in addition to the left-handed one [^3[ p5| . 

First we introduce the model. This model has three families of matter multiplets tpi, <pi, and r]i, which 
are 10, 5*, and 1 dimension representations of SU(5), respectively, ipi contains the quark doublet, the 
charged lepton singlet, and the up- type quark singlet, while 4n the down- type quark singlet and the lepton 
doublet and rji the right-handed neutrino, respectively. The model has 5 and 5* dimension representation 
Higgs multiplets, H and H . H consists of the MSSM Higgs multiplet H2 and a colored Higgs multiplet 
He, and H another MSSM Higgs multiplet Hi and another colored Higgs multiplet He- The GUT gauge 
symmetry is spontaneously broken into the SM one at the GUT scale Mqut — 2 x 10 16 GeV. Above the 
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GUT scale the superpotential W of the matter sector of this model is 

W = -^ u ^ B ^ D H E e A BCDE + V2f dij ^ B ^AHB 

where A, B, ... arc indices of SU(5) and run from 1 to 5. We also introduce the soft SUSY breaking terms 
associated with the GUT multiplets. The relevant part of them is Q 

- £SUSY breaking = (w^)ijV'l^ + ( m 4>)ij<t > i<t > j + ( m l)ijvlvj + mffi h + Vflffih 

+ j^Auy^i^/i + V2A dlj ipi(f>jh + A Uij fji<j>jh + fi,.c. j , (33) 

where ipi, 4n, and f\i are the scalar components of the ipi, (f>i, and rji chiral multiplets, respectively, 
and h and h are the Higgs bosons. In the minimal SUGRA scenario these coefficients are given at the 
gravitational scale as 

m l = m t = m o> 

A Ulj = / Mlj ao, A dij = fd tj a (h A Vi . = f Vlj a . (34) 

At the GUT scale we choose a basis where the up-type quark and the neutrino Yukawa coupling 
matrices are diagonalizcd as 

f = f e^ u i8-- 

J Uij J Ui^ U IJ J 

fdij = (VKM)ikfd k (V^ ) )kj , 

fu ti = U^Sn, (35) 

where f % p i (ip = u, d, v) are the eigenvalues of /,/,.. , respectively, Vkm the Kobayashi-Maskawa matrix 
at the GUT scale, and Vd a unitary matrix which describes the generation mixing in the lepton sector. 
(fii/jiiip — u,v) are phase factors which satisfy 4> Ul + (f> U2 + <p U3 ~ and (p Vl + 4> V2 + 0„ 3 = 0. However 
these phases are completely irrelevant for our below discussion. 

At the GUT scale the Yukawa coupling constants responsible to the down-type quark masses and 
those responsible to the charged lepton masses are supposed to unify as 

(36) 

4 For simplicity we neglect the Yukawa coupling XHT,H and the soft SUSY-breaking parameters associated with it, where 
£ is an adjoint representation Higgs multiplet causing the breaking SU(5)gut - * SU(3) c x SU(2)lX U(l)y. 
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This relation is consistent with the particle spectrum at the low energy only for the third-generation. In 
order to explain the fermion masses of the first- and the second-generations, one has to consider the effect 
of the nonrenormalizable terms also. At that time those terms can be another source of LFV |2^, |2?J , 
but we do not take them into account for simplicity. 

Below the GUT scale we take the basis in which the Yukawa coupling constant matrix responsible for 
charged lepton masses is diagonalized. The basis we take at low energy region is related to that of GUT 
multiplets by the following embedding: 

4>i = {VDijDj,VDijLj}, 

m = {e-^Ni}. (37) 

Then the superpotential W is expanded in terms of the MSSM fields as 

W = f Ul QiU l H 2 + {V^ M ) l0 f d] Q l 'D, J H 1 
+fdiEiLiHi — f Ui VDijNiLjH2 
+fuj{V\<iM)jiEiUjHc — -fui el ^ Ui QiQiHc 
+(V£ M ) ij f dj e- i ^U i D j H c - (V&JijfaQiLjHc 
+f Vi VDijNiDjHc 

+^M UiUi N i N j . (38) 

Here we should notice that the fifth term of the right-hand side of the above equation is no longer 
generation-diagonal. This is nothing but a direct consequence of the GUT unification, that is, one of 
the central goal of the grand unification is to embed the leptons and the quarks into the same multiplet, 
which forces the mixing in the quark sector related to that of the lepton sector. No redefinition of E 



in generation-space can eliminate this mixing, as can be seen from Eq. (38). This mixing causes the 
off-diagonal elements of (mf ) via radiative corrections. 

As for the origin of the observed mixing angle between the left-handed neutrinos a parallel discussion 
to that of the previous section applies. The Majorana mass matrix in Eq. ( |38| ) has an inter-generational 
mixing as 

M ViVj =Ut k M Vk Uty (39) 
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The mass matrix of the left-handed neutrinos (m„) is then 

(m v )ij = Vn ik (m v ) k iV D ij, (40) 

where 

= V Mtk m v k V Mkj- (41) 

Here also m Ui D = fvfV sin / \2, the same notation as that of the previous section. The discussion in the 
previous section shows that the large mixing angle from Vm requires a fine-tuning between the elements 
of M UiVj (Eq. (|l2|)). Therefore as for the large mixing angle between neutrinos it is natural that its origin 
is in the mixing matrix Vd in Eq. (|3^). Here we assume Uij = for simplicity, which means that the 
mixing comes only from Vd- 

Now we evaluate the off-diagonal elements of the slepton mass matrices at the low energy. As stated 
above, both the left- and right-handed slepton's ones have non-negligible off-diagonal elements at one- 
loop level. Assuming m Ue <C OT^ <C m VT we neglect f Vl , and also f Ul and f U2 to obtain approximate 
formulas for the off-diagonal elements of the slepton mass matrices as 



(r4h - -^2/ u 2 3(^KM)3 l (^KM)3 J (3m2+ag)log I ^, (42) 



K)« - -^2 [foVBuVrnj log = + flV£ 2i V D2j log ^ j (3ml + og), (43) 



(44) 



for i 7^ j. These formulas are obtained by a logarithmic approximation from the RGE's (given in 
Appendix |d|). 

Now we study the individual LFV processes. First we concentrate on the r — > fi'y decay. For r — > /i7 
the most important contribution is from diagrams which involves (m|)32 and winos (that is, diagrams of 
Figs. ( |l0| ) (a) and (b)), which is the common feature with the MSSMRN case. These diagrams dominate 
because the large Vi>32 element, suggested by the atmospheric neutrino anomaly, enhances (m?)32 as 

(m|) 32 ~ -^(3mg + al)V^V m2 fl 3 log (45) 

which is the same situation as in the MSSMRN case. The main difference from the MSSMRN case is 
a presence of (rrie)32, but the contribution to the r — > fij is too small at the broad parameter region 
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to be comparable to those from Figs, (jig) (a) and (b), because (mg)^ is suppressed by small (Vkm)32 
p8| . Our result of numerical calculation, Fig. ([ll]), indeed shows that almost the same situation as in 
the MSSMRN case is realized. In the figure we plot the dependence of the branching ratio of t — > /ij 
on the third generation right-handed Majorana mass M„ 3 for tan/3 = 3, 10, and 30. The upper curve 
corresponds to larger tan/3. We take the bino mass as 65GeV, the right-handed selectron mass 160GeV, 
and the tau neutrino mass 0.07eV, as expected from the atmospheric neutrino result. We take ao = 
for simplicity. The figure shows us that the branching ratio of r — > /17 is nearly proportional to the 
square of M Va . At the right-hand side of each curve the Yukawa coupling constant f„ 3 blows up below 
the gravitational scale, so the perturbative treatment is no longer valid in this region. In the region 
near M Vi ~ 10 14 GeV the branching ratio is close to or even beyond the current experimental bound, 
Br(r -> ^7) < 3.0 x 10~ 6 @. At relatively small M„ 3 region (M„ 3 < 4 x 10 12 GeV) the contribution from 
the right-handed slepton mass matrix (the diagrams shown in Figs. (^0|)(c) and (d)) tends to dominate, 
and the curves of the branching ratio show deviation from simple straight lines. 

Next the /i — > process. Since we know Vd32 — C(l) from the atmospheric neutrino result we can 
calculate (m?) 2 3 as 



K>3 - -^Z2fu 3 V^ 2 V D33 (3m 2 + al) log (46) 



g-^/'a VSaaVbsa (3m^ + a%) log ^ 
On the other hand {m\) 3 \ is determined from the GUT symmetry as 



(mf) 31 ~--^/2 3 (Vk M )33(^M)3i(3mg + a g)log^-. (47) 

Then we can definitely calculate the diagram shown in Fig. (|i~2|)(e), which contributes to A^^ defined in 
Eq. (p0|). This contribution is proportional to m T , and it tends to dominate over the other contribution 
to A^^ . We need to know Vd31 an d Vd21 to calculate A^^ ■ However, we can evaluate only the lower 



bound of the branching ratio from the structure given in Eq. ( |21| ) even if we do not know them. The 
contribution from the diagram (e) can be so large that it reaches the present experimental bound at some 
parameter regions |Pq| . 

We can imagine some cases where much larger rate than this lower bound is predicted by the contri- 
bution from A^ e \ One of such cases is that Vd31 is large. In this case, the diagrams in Figs, jl^ ) (a)-(d) 
and (f) are enhanced since (771^)31 and (m|)2i are proportional to VD31 as 

K)3i * -^( 3m o + a o)/" 3 VSasVbsi log (48) 
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K)2i * -^(Smg + ^VSaaVbsilog-jP. (49) 
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While the diagram (f) is enhanced by m T , the contribution is suppressed by (m§)32, which is proportional 
to (Vkm)32- As a result, the diagrams (a)-(d) dominate since Vd32 is large. Also, if /„ 2 and Vd21 are 
non-negligibly large, (m|) 2 i is enhanced as 

(m|) 2 i * -± (4^32^31 log ^ + £ VfaVim log ^) (3m 2 + og), (50) 

and the diagrams (a) and (c) dominate over the other contributions. 

Now we examine these expectations numerically. First we set f U2 to zero and later investigate the 
non-zero f U2 case. We calculate the dependence of the branching ratio of /i — > ej on M U3 and Vdsi in 
the Figs. ( p^| , p| ). In the figures we take V2532 = — 1/V2 and the tau neutrino mass 0.07eV, as suggested 
by the atmospheric neutrino result. We neglect m v here. Other parameters are the bino mass 65GeV, 
the right-handed selectron mass 160GeV, tan/3 = 3 and 30, and the Higgsino mass /i > 0. From the 
Figs. (|l^JTJ) we can see that for Vd31 ~ 10~ 3 the diagram (e) dominates over the others. At relatively 
larger Vdzi region (Vd 31 > 10- 2 - 5 ) the dia grams (a)-(d) become dominant, and the predicted rate is large 
enough to reach the experimental bound for V D31 ~ 10~ 2 if tan/3 = 30(3) and M V3 > 10 13 ' 5 (10 145 )GeV. 

Next let us consider the finite /„ 2 case. We show in Figs. (EM |l6|) the dependence of the /i — > 
branching ratio on Vd31 ancl t ne typical right-handed neutrino Majorana mass Mjy. In the figures the 
parameters that describe the MSW large angle solution are taken as the same as those we used in the 
MSSMRN case, and other parameters are taken as the same as in Figs. (13, 14|). We assume in the figure 
the universality of the right-handed Majorana mass, that is, M Vx = M V2 = M V3 (= Mn). We can see from 
the figures that the enhancement due to the MSW large angle solution is so large that the dependence 
of the rate on Vd31 is small, and that the almost same situation as in the MSSMRN is reproduced here 
again, as can be seen when comparing these figures with Fig. (§). Since the dominant contribution is 
determined by the second term of Eq. (p0|), the rate is almost determined from the value of M„ 2 because 
/ 2 2 cx M„ 2 for fixed m V)i . Hence we can conclude from the Figs. ( |I1^ , |l6| ) that for tan/3 = 30(3) the 
excluded region of M„ 2 extends to 9 x 10 12 (1 x 10 14 )GeV, at least for the parameters chosen in our 
calculation. 

We also calculated for the parameters suggested from the MSW small angle solution and the 'just so' 
solution. For the MSW small angle solution the difference from the f V2 = case is at most the factor 2 
enhancement, and for the 'just so' solution no difference from the /„ 2 = case can be seen. 

Finally we comment on the b — > sj process. This model has a characteristic feature that the right- 
handed neutrinos couple to the right-handed down-type quarks with the large mixing through the ninth 
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term of Eq. (pq ) . This coupling induces the off-diagonal elements in the right-handed down-type squark 
soft mass matrix via the radiative corrections, which are as large as 

(mf )ij * -^f*.VB kj V Dki (tor% + a 2 ) log ^L. (51) 

This fact causes us an expectation that the b — > sj rate may become larger by the effect of the diagram 
which contains (m~ )23, by an exact analogy to that of the lepton sector. We examined whether this is 
true or not, and concluded that a tiny enhancement indeed occurs but it is too small to experimentally 
distinguishable. This is because the heavy gluino suppresses the contribution. 

4 Conclusion 

In this article taking the solar and the atmospheric neutrino experiment results into account we in- 
vestigated the lepton flavor violating decay processes, such as /i — ► or t — * /ij in the MSSM with 
the right-handed neutrinos (MSSMRN) and in the SU(5) SUSY GUT with the right-handed neutrinos 
(SU(5)RN). 

In the MSSMRN we first studied the branching ratio of r — > fjrf. It gets larger for larger tan/3 and 
almost proportional to M„ for fixed tyi Vt . A large rate naturally follows from the large mixing between 
v T and i^j, suggested by the atmospheric neutrino result. If tan ,9=30(10) the rate reaches the current 
experimental bound for M V3 ~ 2x 10 14 (6 x 10 14 )GeV. We investigated the fi — > ej rate under three kinds 
of the solar neutrino solutions, the MSW large and small angle solutions and the 'just so' solution. We 
argued that the [i — ► ej depends on /„ 2 , and that especially the large Vd21 and the large f V2 , which the 
MSW large angle solution suggests, naturally results in such a large rate that the predicted rate is beyond 
the current experimental bound if M„ 2 is larger than 8 x 10 12 (8 x 10 13 )GeV for tan/3 = 30(3). We also 
investigated the dependence of \i — > on Vd31- For tan/3 = 30 and Vd3i = 10~ 2 , M U3 ^ 3 x 10 13 GeV is 
excluded for our input parameters. 

In the SU(5)RN we calculated the t — » /ry rate. Here also the large mixing between v T and leads 
to a large rate, which is almost the same situation as in the MSSMRN. For fi — > e-f in this model, we 
can predict the lower bound of the branching ratio of it. This lower bound is calculable from the value 
of m„ T and the large mixing angle between v^v^, suggested from the atmospheric neutrino result, and 
it turns out to be within accessible region by near future experiments. They are supposed to probe the 
fi — > ej branching ratio to 10~ 14 level [p|, and then the region M V3 > 10 13 (10 12 )GeV can be probed 
for tan/3 = 3(30) and m Vr = 0.07eV. We considered the relation between fi — * ej and three kinds of the 
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solar neutrino solutions. The large rate naturally follows from the MSW large angle solution, similarly 
to the MSSMRN case, while in the MSW small angle solution and in the 'just so' solution there is only 
a little difference from the f V2 = case. 

If the LFV processes are discovered or the experimental bounds are improved by near future exper- 
iments, the interesting insight on the lepton sector will be obtained. The best effort to implement it is 
strongly desired. 
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A Definitions and Conventions in the MSSM 



In this appendix we collect our notations of the MSSM used in our article. The superpotential of the 
MSSM Wmssm is defined as 

V^mssm = fe^HtEiLj + f d 1 1. (IP, + /„ HSir, + nHxHi. (55) 

We use i and j as the generation indices running from 1 to 3. E, _D, and U are the superfields associated 
with the right-handed electron e^j, the right-handed down- type quark cIr, and the right-handed up- type 
quark ur, respectively. Q and L are ones associated with the quark doublet and the lepton doublet 
II defined by 

9l=( ^V^=f:. ), (56) 



,4/' V e L . 
and Hi and H 2 the Higgs doublets whose SU(2)^ components we denote as 

Hi \ „. _ ( at 



h ^{h- J^ 2 = U! )■ (57) 

The scalar components of H® and H% develop the vacuum expectation values (vev's) as 

<B?>s^, (H° 2 ) = ^, (58) 
which satisfy v\ + v\ = v 2 with v ~ 246 GeV. We define the ratio of these vev's as tan j3, 

tan/3=— . (59) 

Vl 

The Yukawa couplings are given by the fermion masses and the Kobayashi-Maskawa matrix as 



i— nip- „ 
v cos p 

v cos p 



/«« = ^-^(Vk M )«. (60) 
u sin p 



Supersymmetric couplings and masses of chiral multiplets are given as 

-C= d 2 9 W + h.c. (52) 



Our convention is unusual in order to keep the chargino mass matrix in accordance with the Haber-Kane convention |29| . 
6 We implicitly assume the contraction convention over SU(2)^ doublet indices (a, b, .. = 1, 2) of two doublets A and B 



(53) 



AB = e ab A a B b , (54) 



where e a (, is an antisymmetric tensor with e\2 = — €21 = 1. 
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We also introduce the soft SUSY breaking terms in the Lagrangian, 

- C susy breaking = {rn\)ijV Li lLj + (mj) ij e* Ri e Rj 

+ { m %)i]^ L SL 3 + {ml)ijU Ri u Rj + {mjjijd^dRj 
+m 2 hl h\hi + m\ 2 h\h2 

+ B h hih 2 + h.c.) 
+ (^M 1 B L B L + X -M 2 W a L W a L + ^M 3 g a L 9l + h.c). (61) 

Here the first seven terms are the soft SUSY breaking masses for the doublet-slepton I l , the right-handed 
charged slepton e.R, the doublct-squark qx, the right-handed up-type (down-type) squark ur (c?k), and 
the Higgs bosons. B, W 7 and g stand for bino, wino, and gluino respectively, and the superscript a is the 
gauge group index for each corresponding gauge group. 

Now we discuss the slepton mass matrices. Let e^i and e.Ri be the superpartners of the left-handed 
electron e Li and the right-handed electron e Ri , respectively. Then the slepton mass matrix (m 2 )ij is 

- £ -(4.4)^(«-)-(4.4)( ( S& ttO(S)' <«> 

Here m\ and m R are 3x3 hermitian matrices and m 2 LR is a 3 x 3 matrix. They are given as 

( m l)ij = ( m |)ij + m l^i] + m%5 lj cos2/3(-i + sin 2 9 W ), (63) 
(m|j)y = (mf)jj + m^Si-j - m 2 z 6ij cos 2/3 sin 2 9 W , (64) 
( TO lfl)y — A\ j v cos (3 / V2 — m ei [i5ij tan (3. (65) 

As for the neutrinos we should notice that there is no right-handed sneutrino in the MSSM. Let i>i be 
the superpartner of the left-handed neutrino Vi. The mass matrix of sneutrino (m 2 )ij is 



2 1 

l l>ij + 2 nu Z u iO 



(ml)ij = {m 2 j)ij + \m%5ij cos 2/3. (66) 



Next we discuss the mass matrix of gauginos. First we consider chargino mass matrix Mc- It is a 
2x2 matrix that appears in the chargino mass terms, 
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Mc is diagonalizcd by 2 x 2 real orthogonal matrices Ol and Or as 



R M c Ol=dmg(M--,M.-). 



Define the mass eigenstates Xal(r) = 1, 2) by 



x l L ) = o L 

X 2 L ' 



H 



1L 



Xl - R ) = Or 
X 2 R J 



W 
H. 



R 



2R 



Then 



(A =1,2) 

Finally we consider neutralinos. The mass matrix of the neutralino sector is given by 



Xa = Xal + Xar 
forms a Dirac fermion with mass M~- . 



-c m = \ (b l w q l h° l h 2L ) M n 



( B L \ 

wl 

H° 1L 



+ h.c. 



where 



/ 



M N = 



Mi 






Mo 



— mz sin 9w cos (3 mz cos 9w cos (3 
\ m z sin 0\y sin [5 — mz cos 9w sin (3 

The diagonalization is done by a real orthogonal matrix On 



—mz sin 6 'w cos j3 mz sin Ow sin [3 \ 
cos cos /? — TOzcos^vKsin/3 
-/i 








N M N N = diag(M^o, ...,M^o). 



The mass eigenstates are given by 



Xal ~ {On)abX° bl 



(A,B = !,-■-, A) 



where 



We have thus Majorana spinors 



x al - (Bl,WI,H^ l ,H% l ). 



Xa = Xal + Xar, [A = 1,-", 4) 



with mass M^o 



(68) 



(69) 



(70) 



(71) 



(72) 



(73) 



(74) 



(75) 



(76) 



Now we give the interaction Lagrangian of lepton-slepton-chargino (-neutralino) in a basis of the 
slepton weak-eigenstate and the chargino (neutralino) mass-eigenstate. By writing the interactions in 
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this basis we can get transparent view in the discussion of the multimass insertion technique discussed 
in Appendices B and C. The Lagrangian is 

-Ant = vlrAiC^PR + C^P^e, 

+4 ¥ A {N^Pr + N$£>P L ) e t + h.c, (77) 



where the coefficients are 



C A L f = 9 2(0 R ) A1 , 
m V2m ez 

C LR — ^{^L) A 2, 

V COS p 

n ll = ^=[-(O N ) A2 -(O N ) A1 t & n0 w }, 
V2 



v cos p 

= V2g 2 (0 N ) A1 tnn9 w . (78) 

B Multimass Insertion Technique 

Mass insertion technique is useful to understand lepton flavor violating processes in the MSSM since only 
off-diagonal elements of the slepton mass matrices are sources of the flavor violation. In this appendix 
we introduce the multimass insertion technique with the nondegenerate masses. The multimass inserted 
diagrams may not be necessarily suppressed than the single-mass inserted ones when the slepton mass 
matrix has several small flavor-violating elements. Also, though in the previous mass insertion formulas 
the degeneracy of the slepton masses is sometimes assumed, such a degeneracy is not necessarily main- 
tained at low energy even if the universal scalar mass hypothesis is assumed in the higher energy scale. 
Our rule to derive the mass-inserted amplitudes is very simple. We can derive them by taking the finite 
difference on amplitudes which have no mass-insertion. In this section we derive general formulas of 
e\ — ► e^7 keeping in mind that we apply them to more realistic cases of /i + e + j and r + — > /i + 7 in 
the next section. 

We refer to the slepton mass matrices as (rn~), 

-£= E i m W\h (79) 
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where 

h = (eL,ML,f L ,e fl ,,/2 K ,f fl ) T (for lj = e t (i = 1 - 6)), , g0 > 
Vi = {v ei v^v r Y (for li = Vi (i = l-3)). 

The explicit forms of these matrices are given in Appendix Here we assume that all the off-diagonal 

elements of the slepton mass matrices, both flavor violating and conserving ones, are much smaller than 

the diagonal elements, ((fn~)u ^> (rn~)jk (j ^ k)). 

In the mass-insertion technique the internal slepton lines are classified to two types at one loop level. 

First type is a slepton line on which the momentum of slepton is not changed as 



h 



>N 



(mf) h i 2 {mf)i 2h (mf) hh (mf)i N _ lh 



^ir^^^^KWa^rK)^ • • • {m l ) lN _ xlNW - w , (81) 



where we call the diagonal components of the slepton mass matrices, (rrt~)u, as rn~ . We would like to 
consider the — > e^j process, and so in the above figure we mean that an anti-slepton is going from 
left to right with a momentum k. The product of propagators in above equation is referred to as i*jy, 



Fjvfm- ,m~ ,m~ ) = "_ ? — — ■ (82) 

h h In' k 2 -m 2 k 2 -mf k 2 -rn 2 



I 1 1 

i~ k 2 — m~ k 2 — ?.„, 

The second type is a slepton line where the momentum is changed, by an emission of a photon with an 
outgoing momentum p, as 




h h h In 
• • • 

(m 2 ) hh (m?) /2 ; 3 (m 2 ) hh (M 2 )i N _ 1 i N 



p + k 
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= (k + P y-ml k 2 -ml K)w 2 ^3=^K)«3 • • • K) iN _ liN ^=r- 

+ (k + p) 2 - ml {m l )hh (k + p) 2 - ™l fc2 - K ^ )/2i3 ' ' ' {m t )lN - llN k*-ml 
+ ■■■ 

(k+p) 2 —ml^ ^ llt2 (k + p) 2 — rf%l^ '^ 2 ' 3 ^ ^ 1n ~ i1n (k + p) 2 — mf N k 2 — mj N 

(83) 

The sum of product of propagators in this equation is referred to as Gat, 

,_2 -2 -2 > 1 11 1 



Gjv(mr ,m.- ,"-,m,- ) = 



^'"V " (k+p) 2 -m 2 k 2 -m? k 2 -mf k 2 -mf 



h h h In 

I 111 

(A; + p) 2 - ml (k + p) 2 - m 2 ^ k 2 - m~ k 2 - m 2 ^ 
+ ■■■ 

II 11 

+ 



(k + p) 2 - m~ (k + p) 2 - m~ (k + p) 2 — m~ k 2 — to- 



(84) 



These functions Fat and Gat can be given as a finite difference of Fn-i and Gn-i, 



F N (m Ii ,m !2 ,---,m fN ) = D^at.^to^to^, • • • , m^); mj]{m- h , to^), (85) 

G A t(to? i ,to? 2 ,---,to? jv ) = D[G N - 1 (m 2 I ,m? 3 ,---,m 2 J;m 2 l ](ml,m'f 2 ), (86) 

where 

D[f(x);x](x 1 ,x 2 ) = — L_ [f{xi) - f(x2)\ . (87) 

Then, by taking finite difference in sequence, each scalar line can be represented as a linear combination 
of flavor conserving scalar lines, F\ or Gi, 

F N (m 2 i ,m 2 2 ,---,m 2 N ) = D N - 1 [F 1 (m 2 y,m 2 ](m 2 i ,m 2 2 ,- ■ ■ ,mf N ), 
G N (ml,ml,---,m 2 N ) = D N - 1 [G 1 (m 2 );m 2 }(ml 1 ml,---,m 2 iN ), (88) 

where 



D N ~ 1 [f(x); x](xi, X2, ••• , Xn) = Elll^— )f( x i)- 



(89) 



Due to this fact, we can get amplitudes with any masses inserted from the corresponding flavor-conserving 
diagrams. 
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Before we derive amplitudes of el — > e^7 (i > j), we introduce a mass function iff as 

/" r1 4 k N 1 

^K.mi.^m^) = y^]^n ] _L_ f . (90) 
This function ijjf can be reduced to IJjjLi or T^f -1 by following rules, 

In (ml, rn%, ■ ■ ■ ,m%) = Dllff^m 2 , ml, ■ ■ ■ ,m N );m 2 )(ml,ml), (91) 

lff(ml,m%,---,m N ) = m N lff~ x (ml,m%,- ■ ■, m%) + lffZi(m\ ,m\, ■ ■ ■ ,m\_ l ), (92) 

and then, all iff for N > 2 + M can be derived from I®, 

/iV 2 ) = (l-log^)m 2 , (93) 

where A is a renormalization point. Also, the signs of Iff and D L [lff] are definitely determined as 

(-l) N + M lff(ml,ml,---,m N )>Q, (94) 
(-l) L + M+N D L [lff( m 2 ,m 2 , ■ ■ ■ , to 2 , m\, ■ ■ ■ , m 2 p ); m 2 ](M 2 , M% , ■ ■ ■ ,M| +1 ) > 0. (95) 

iV-P 

Then, we can discuss about relative signs between several diagrams definitely by using iff and the finite 
differences. Furthermore, since the mass dimension of iff is (4 + 2M - 2N), we can derive a following 
relation, 

N 

(2 + M-N)lff(m 2 ,ml---,m 2 N ) = £ ro, 2 I# +1 (m?, ■ ■ ■ , to 2 , to 2 , • • • , m%) (96) 
for iV > 2 + M, since 

^{x- 2 - M+N lff{xm 2 ,xm 2 ,---,xm 2 N )} = 0. (97) 

Due to Eqs. ( ^I| , p2| , p6[ ) there are some ways to represent one function, for example, 

/ 2 (to 2 ,to 2 ,to 2 ,A7 2 ,M 2 ) = -3m 2 7 5 1 (m 2 , m 2 , to 2 , m 2 , M 2 ). (98) 

We will derive amplitudes of e+ — > e+7 (z > j) by the mass-insertion technique. The amplitude is 
generally written as 

T = e e a * Vi(p) ia am (A<£ j) P L + P R ) Vj (p - q). (99) 
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Here, e is the electric charge, e* the photon polarization vector, Vi and Vj the wave functions for the 
external leptons. Assignment of momenta of the external fields is shown in Fig. (^). Since above term 
is violating the lepton chirality, it is convenient to decompose and A^ as 

A™ = A^\ cl +A^\ c2 + A^\ nl + A^\ n2 + A^U 

= 4^1 + 4^2 + AM \nl + A™ \n2 + 4 j) l«3- (100) 

A^fji\ n i tC i come from diagrams in which the lepton chirality is flipped on the external lines (Figs. (a) 
and (b)). jjj| n 2,c2 are contributions from the diagrams in which the chirality is flipped on a vertex of 
lepton-slepton-neutralino (-chargino) [Figs. (|l8|)(c) and (d)]. A^ R \ n 3 are contributions from those with 
the chirality flip on the internal slepton line (Fig. (|l8|)(e)). Subscripts c and n represent that they are 
from chargino and neutralino diagrams, respectively. 

The neutralino contribution to A^f 1 \ n \ , derived from a diagram where the chirality of lepton is flipped 
in the external lepton line, is given by 



N=l h,-,lN-l 

-2 



A 



m i 

R l«l 



x [I£ (mi , m£, mf , M~o , M\ o ) , mg] (m 2 ^ , , • • • , m% N ) 
(Z = i,In = j), 

UfXi) , (ioi) 

V / L<^R,l„=i+3,lN = i+3 



where m ei is the i-th generation charged lepton mass, and the explicit form of /| is 

J 5 2 (mim?,m?,M? ,M?o) = n „ 1 a (1 - 6a; + 3a; 2 + 2a; 3 - 6a; 2 log i) (102) 
A A to| 2(1 — a;) 4 

with x — M? /fn^. is a coupling constant between slepton and neutralino, and our definition of 

couplings of slepton to neutralino (or chargino) is 

+4 ^(N^Pr + N^Pl) e, 

+4 &(<rPr + KlPl) * + h.c. (103) 

Here, , and N^j^ correspond to the (lepton-chirality conserving) gaugino interactions, and 

C'lr 1 ' ^rli an( l ^lr 1 are * ne (lepton-chirality violating) Higgsino interactions. These coupling con- 
stants are represented by the MSSM parameters in Appendix 
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Sign of the term with N off-diagonal inserted- masses in A^\ n i is definitely determined by signs of 
coupling constants of slepton to neutralino and the inserted masses, and is independent of the diagonal 
slepton masses and neutralino masses since 

(-l^+WlJl^^TOf.MloA^,),^]^,^,..-,^) > 0. (104) 

A contribution from a neutralino diagram where the lepton chirality is flipped in the vertex of slepton- 
lepton-neutralino is given as 

^ oo 

^ u = -2fctf m k N ™' N £P E E (^^Wixh-KW-xi* 

^ ' N=\ h,-,l N -i 

x D N [I\ {ml , ml , M|o , M?„ ) ; m? ] (m? q , m? j , • • • , m? !jv ) 
(Zo = z,Zjv = j), 

UPU) ■ (105) 



^i? I "2 



The explicit form of J| in Eq. (105) is 

Jl (m? , ml , M? , M? ) = - ^ n j ^ 3 ( 1 - x 2 + 2x log x) (106) 

/Tig ^1 

with x = M| /mf, and the sign of D^I^] is (-1)^+* from Eq. @. 

A contribution from a neutralino diagram where the lepton chirality is flipped in the internal slepton 
line is given as 

x D N [I\ {ml , ml M|o , M|o ) ; mf] (mf io , mf (i , • • • , m 2 ^ ) 



A m 



R \n3 



{l a = i + 3,l N =j), 

UPlm) ■ (107) 



A contribution from a chargino diagram where the lepton chirality is flipped in the external lepton 
line is given as 



i 00 
A ^ J) ' C1 = 6f47r) 2 TOe ' C ll* c ll E E (^I)'oii(^)/i/2---(wI); JV _ 1 ; JV 

x [I 5 2 (m? , mf , M?_ , M?_ , M?_ ) ; mf] (mf , o , mf ;i , • • • , mf , ^ ) 
(/ =i,l N = j), 

A^\ cl = 0{m ej ). (108) 
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The explicit form of if in Eq. (108) is 



II (ml ,mlM*, M|_ , M|_ ) = - ^ — 1 (2 + 3x - 6x 2 + x 3 + 6x log x) (109) 
* * x m~ 2(1 — xj 4 

with x = M 2 ._/ml, and the sign of D N [lg] is (-l) N+1 . 

A contribution from a chargino diagram where the lepton chirality is flipped in the vertex of lepton- 

sneutrino-chargino is given as 

^ ' N=l h,-,l N - t 

X A X a x A h uiK "to' «tl' ' In 
(la = i,l N = j), 

4 J) |c2 - 0(m ej ). (110) 



The explicit form of 1\ in Eq. (110) is 



I\ (ml ,M 2 ,M 2 ,M 2 ) = ^ — 1 .„ (3 - Ax + x 2 + 2 log x) (111) 



ml 2(1 -xf 



with x = M 2 _/ml, and the sign of D N is (-l) Ar+1 . 

C Application of Mass Insertion Formulas to r + — > /i + 7 and fi + — > 

e + 7 

In this section we apply the formulas derived in the previous section to t + — ► fi + j and fi + — > e + 7 
processes. We neglect m^ (m e ) in the calculation of r + — > /i + 7 — * e + 7). First we consider r + — * /i + 7. 
The whole contribution to t + — > /i + 7 can be written as 

4 TM) = 4 M W4 M W4 M) i*i+4 T/1 W4 T|0 u 

4"> = 4^ ) | c1 +4^ ) | c2 +4^ ) i„ 1 +4™ ) i»2 + 4 ai) U3, (H2) 

in the same notation as the previous section. In many models the dominant contributions to t — ► fj.'y 
are from the diagrams with single insertion of (m-)^, (m|)32, A 32 , or A 23 . Among these lepton flavor 



violating coupling constants only (w~)32 and (m 2 )s2 are important for tan/3^ 1. We here show explicit 

• > 



expressions of the diagrams with single off-diagonal slepton mass matrix element insertion for tan/3^ 1. 
From the formula we derived in the previous Appendix, the contributions from diagrams with the lepton 
chirality flipped in the external line (Figs. (|l8|)(a) and (b)) are given as 

/i( t m)i ^ a2 (n \2 
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4 TM) i-i 



1 K I™ 1 



(m|) 32 


( 


_m| M -ro? r _ 











1 0.2 

mr^-r {{0 N ) A 2 + (O n )ai t&nOw) 

1Z 47T 



( m |)32 
9 9 
M-L TL 

^3 4^ (OAr) ^ 



(™p32 



— fnl(xAjl L ) —fnl{xAf L ) 



m% 



— fnl{XAiX R ) —fnl(xAf R ) 



(113) 
(114) 

(115) 
(116) 



where x A j — M~. A /m~ with I = , v T , [il, tl, M-R, and tr and = x A and \a- The coefficient 
functions / c i(x) and f n i{x) are given as 

1 



/cl(z) = 



: (2 + 3a; - 6x 2 + x A + 6x log x), 



2(1 -x) 4 

fni(x) = 1 (1 - 6x + 3x 2 + 2x 3 - 6x 2 logx), 



(117) 
(118) 



2(1 - a;) 4 1 

which are positive definite and monotonically decreasing. Following coefficient functions are also defined 
to be positive definite and monotonically decreasing. 

The diagrams where the lepton chirality is flipped in the vertices give the following contributions, 



L Ic2 



A R \c2 



A R \n2 



a 2 M 



47r ^/2mw cos (i 

2 



(O fl ) A i(0 L ) A2 



(m|) 32 



1 a 2 

— TOt--- 



M x3 



4 47r to^ cos 6>vk cos (3 

,2 



^-/c2(xap,J - ^-/ c2 (a;ApJ \ , 



(O n ) A3 ((O n )a2 + [O n )a\ t&n6 w ) 



(to|) 32 



— fn2{XAjl L ) \-fn2{xAf L ) \ , 



TO- 



m% 



(119) 
(120) 

(121) 



1 ay 



2 47r to 2 sin 0w cos /3 

32 



Al 



(™p3 



2 2 

to ~ — m% 

HR TR 



— fn2{XAp, R ) —.fn2{xAf„) 



TO ~ 
HR 



where the functions / c2 (x) and f n 2(x) arc defined as 

1 



fc2(x) 
fn2(x) 



(3-4x + x 2 + 2 log a;), 



2(1 -a;) 3 

:(1 — a; 2 + 2a; log x). 



(1-x) 3 ' 



(122) 

(123) 
(124) 



2<S 



Finally the contribution from the diagrams in which the lepton chirality is flipped on the internal 



slepton lines are 



^l M) |«3 = -7}-^(°n)ai ((O n )a2 + {0 N )Aita.n9 w )tane w (m 2 LR ) 33 (m 2 L ) 3 2M^o A 



fn2 {xAf E ) 



1 1 f 

9 9 9 9 9~ 

777? r - m*. tt7? r - 

1 1 1 , V 

2 2 2 2 2 J" 2 \ X At l ) 

m% m% — m% m% — m- 



-fn2(x A fi L )}, (125) 



2 2 2 2 2 

TOr — 777- 777^ — 777| 

Ml. ML r -R ML t l 

^ M) |n3 = -\^{°n)ai {{O n )a2 + (0 N )Ait&n9 w )ta,n8 w (m 2 LR ) 33 (ml) 32 M S( a A 
1 1 1 



2 2 2 2 2 

777^ 777j — 777;= 777- — 777^ 
TL TL T R Tl fl R 



fn2(xAf L ) 



2 2 2 2 2 Jn2\XAf R ) 

m f R m f R ~ m i L m i R - m i R 



2 2 2 2 2 

TO- 777- — 7Tl% 777- — 777- 
MB MB r t MR TR 



fn2(x AilR ) . (126) 



Next we present the mass insertion formula for fi + — ► e + 7. Here we assume mass degeneracy between 
sclcctron and smuon, then 



"II' 



The contributions from diagrams with the lepton chirality flipped in the external line (Figs, 
and (b)) are given as 

„(pe), la 2 



i)(a) 



Id 



.4 



.4 



(Me) I 



2 

AX 



(raf ) 2 i + 



(m 2 f ) 23 (m 2 f ) 3 i 



f 



rPci(^Ai?) 



(m|) 23 (m,|)3i 

2 \2 



(777? - 777? ) 



-fcl(XAO T 



0(m e ), 



1 a 2 



~ m M7^T~ ((Oat)a2 + (Otv)ai tanflvt/)" 

iZ 47T 



(m?) 2 i + 



(m?) 2 3(m?) 3 i 



777? — 777? 



l( X Alr) 



m 4 ym\x A j L 



(127) 
(128) 
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(m|) 23 (m|) 3 i 



(to? — TO? ) 2 



3 m ^(0w)Ai 



"4 



TL 



(129) 



("»I)21 + 



(m|) 23 (m|)3i 



4 

TO- 



(to?) 23 (to?) 3 i 
(to? h -to?J 2 



TO? 



fnl{xAf B 



(130) 



with 



5ci(z) ee / cl (x) +x/^(x), 

9ni(x) ee +x/; i (x). (131) 

Here primes on f c i, n i{x) mean differentiation about x. 

The diagrams where the lepton chirality is flipped in the vertices give the following contributions, 

a 2 M-- 
1 47T y/2mw cos P R A1 L A2 
(TO?) 23 (m|) 3 i 



-TO, 



(m|) 2 i 



4) 



2 

TO- - 
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- m~ 






(™|)31 




TO?J 2 





—g C 2(x A v) 
mr 



-2-fc2{XAO T ) 



^ e, | c2 = 0(m e ), 

4" 6) |n2 - 7 -(0 N ) A 3((0 N )A2 + (0 N )Alt & n6 

4 An TOZ cos ffw cos p 

(m|) 23 (m|) 3 i 



(132) 
(133) 
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«Mor ^(0 W )A 3 (07v)Al 

2 An mz sin 6V cos p 



(to?) 2 i + 



(to?) 23 (to?) 31 
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to~ — m| 



—T9n2{x A l) 
Ir 



(w1) 23 (to?) 3 i 

(TO? -TO?J 2 



— 2-/n2(arAf H ) f • 

m TR I 



(135) 



Here the functions g c2 {x) and g n2 (a;) are defined, similarly to g c i,ni(x), as 

g c2 (x) ee f c2 (x) + xf' c2 (x), 
g n2 (x) ee ,f n2 {x)+x,f' n2 {x). 



(136) 
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The contributions from Fig. (|18|)(e) in which the lepton chirality is flipped in the slepton lines are 
following, 

4^ e) |„ 3 = —Afgo [O N ) A i{(O N )A2 + (O N )Aitaii0w)tsm6w 



x <J -(m s )23(m LR ) 3 3(m L ) 3 i 

( 1 1 1 1 



> > > » > > 9 fn^i^Alu) 

mi mi — to| mi — mi mi — mi R 

Ir Ir t R Ir t l Ir II 

1 1 1 1 , _ , 

9 9 9 2 2 2 2 Jn2\p^Afi, ) 



9 9 9 9 9 9 9 

m| mi — to- mi — mi mi — m~ 

TR Tr I Tr T L Tr I 



I 1 1 1 , _ A 

' 2 9 92 99 2 Jn2\X A] T ) I 

to- to- — to- — mi mi — to- l ) 

II II Ir II tr II tl / 

+ (m| H ) 22 (TO|) 21 

/ 1 1 1 1 N 

x 2~ 7 2 2 \9 J n2 \ x Ain ) ~i 4 2 2~ 5«2 (^ 4f r ) 
\ to- to- — To- ) to- mi — mi ™ L 

\ Ir V *r *z/ (i Ir / 

+ (™lfl)22(m|) 23 (m|) 3 i 

xf- * a 1 av> a 1 2 /„ 2 (^) 

\ to? (to? -to?) 2 to? -to? v 



1 1 1 



2 _ m 2^ 2 ^2 _ m 2 JnZ\^AT L 



fn2(XAf 



TL V TL l L ' T L l R 



mf TO? — TO? TO? — TO? " ! ' 
II II Ir II l 



9n2(x Alr )}}, (137) 



( 3 ) 1 ^-2 

|n3 = 2 4^ M x / ,( O w)^i((°A r )^2 + (OAr)Aitan%)tan6' w 
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+ (-L)2 2 (mI) 23 (™I) 31 

* 2~~ 7 2 2 \9 2 2 f n 2\ X Alr ) 

\ II y II W II Tr 



1 1 1 



■fn2{xAf a ) 
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H — 1 T~~ 2 2 — " —9n2(x A r)) \ . (138) 
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D Renormalization Group Equations 

In this Appendix we show the one-loop level renormalization group equations (RGE's) relevant to our 
discussion in the text. In the equations below we use a shorthand notation, that is, for matrices in 
generation space A and B we define {A, B} by 

{A, B}ij = ^(A ik B kj + B lk A kj ). (139) 

k 

D.l RGE's for the MSSMRN Model 

First we devote ourselves to the MSSMRN model. The RGE's for the gauge coupling constants and the 
gaugino masses are unchanged from the MSSM since the right-handed neutrinos are singlet under the 
Standard Model gauge group. First we list the RGE's of the Yukawa coupling constants. 

leTrV^/e,, = {-|9?-3.g 2 2 + 3Tr(/ d /]) + Tr(/ e /t)| /ei3 

+Xfeflfeh+(fe.ftfuh, (140) 

leTr 2 ,^/^ = | - \g\ - ig 2 2 + 3Tr(/„/t) + Tr^/t) } f vtj 

+Z{UftU)ij + {Uftfe)ij. (141) 

Next RGE's of the soft masses and A parameters. 

IGttV^K)^ = - (jjffi l M i| 2 + Q 9l |M 2 | 2 ) <% - |g 2 ^- 

+2 {ftmlU + mlJtU + AlAe).. 

+2 (ftmlU + m 2 h2 flU + AlA v ).. , (142) 
16^^(777?),, = _| 5 2 |Mi|2j .. + 6^ 

+2{mlf e ft} ij 
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+4 (/ e m|/J + 4/e/I + A e At) . . , (143) 
^ 2 ^(ml)ij = 2{ml,f v fl}ij 

+4 + mL/,/J + A v Al).. , (144) 

l6ir 2 ^A eij = {- 9 -gl-3gt + 3Tr(flf d ) + Tr(pJe)}A eij 

+2 j-^Mi - 3g*M 2 + 3Tr(flA d ) + Tr(/ e U e ) J / ei . 

+4(/ e /t^e)< j + 5(A e /t/ e )^ 

+2(/ e /t^).. + (Ae/ t /[/) .. 5 ( i45) 
\^^A Vii = |-^2_3 5 2 + 3Tr ( / t /u)+Tr(/ t /i/) |^ i . 

+2 j-j^M! - 3g 2 2 M 2 + 3Tr(/t J 4 u ) + Tr(/t^)| /„„ 
+2{f v ftA e ) ij + {A v flf e ) ij , (146) 



where 



S = Tr(m£ + m 2 - - 2m£ - ro£ + ro£) - + m^ 2 . (147) 

Here g\ is the U(l) gauge coupling constant in the GUT convention, which is related to the U(l)y gauge 
coupling constant gy by g Y — §g 2 - 

D.2 RGE's for the SU(5)RN Model 

Next we list the RGE's for the SU(5)RN model. The superpotential of the matter sector is 
W = ^f u ^f B ^f D H E e AB CDE + V2f d ^f B ^ A H B 

and the soft SUSY breaking terms are 

- £susy breaking = (m^k^l^ + {m%)ij<i>\<j>j + (m 2 , ) ij t?J ??j + m\h)h + rrffih 
+ ^ A uij4>ii>jh + V2A dij ^jh + A Vij r)i4>jh + /i.e. j 

+^M 5 X 5L X 5L + h.c. (148) 

We denote the SU(5)gut gauge coupling constant as g$, the SU(5)gut gaugino A5, and its soft Majorana 
mass M 5 . We neglect a coupling \HY>H, as stated in the text. 
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First is the RGE's for the dimensionless coupling constants. 
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Next is for the soft masses. 
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(a) (b) 

Figure 1: The Feynman diagrams which give dominant contributions to t + — > fi + j when tan (3^ 1 and 
the off-diagonal elements of the right-handed slepton mass matrix are negligible, as in the MSSM with 
the right-handed neutrinos. In the diagrams, {m 2 ^)z2 is the (3,2) element of the left-handed slepton soft 
mass matrix. tl(r) and £il(r) are the left-handed (right-handed) stau and smuon, respectively, and v T 
and the tau sneutrino and the mu sneutrino. Hi and H2 are Higgsino, W wino. The symbol fi is the 
Higgsino mass. The arrows represent the chirality. 
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Figure 2: Dependence of the branching ratio of t — > /Z7 on the left-handed selectron mass mg^ in the 
MSSM with the right-handed neutrinos. We take the tau neutrino mass 0.07eV and Vd32 = —0.71, which 
are suggested by the atmospheric neutrino result. M„ 3 is fixed at ~ 10 14 GeV by imposing a condition 
fu 3 = fv 3 at the gravitational scale. The dotted line shown in the figure is the present experimental bound. 
We set the wino mass M2 130GeV, and the Higgsino mass parameter /i positive. The mu neutrino mass 
is neglected. We take tan/3 = 3, 10, and 30. The larger tan/3 corresponds to the upper curve. 
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Figure 3: Dependence of the branching ratio of r — > fxj on the third-generation right-handed neutrino 
Majorana mass M„ 3 in the MSSM with the right-handed neutrinos. The input parameters are the same 
as those of Fig. (2) except that in this figure we take mg L = 170GeV and that we do not impose the 
condition f U3 = /„ 3 but treat M„ 3 as an independent variable. The dotted line shown in the figure is the 
present experimental bound. Here also the larger tan/3 corresponds to the upper curve. 
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Figure 4: Candidates of the Feynman diagrams which give dominant contributions to fi + — ► e + 7 when 
tan (3 ^ 1 and the off-diagonal elements of the right-handed soft mass matrix are negligible, as in the 
MSSM with the right-handed neutrinos. In the diagrams, (m%)ij is the element of the left-handed 
slepton mass matrix. &l(r) is the left-handed (right-handed) selectron and v e is the electron sneutrino. 
Other symbols are the same as those in Fig. (1). 
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ji^ey in the MSSMRN with the MSW large angle solution 




M V2 (GeV) 

Figure 5: Dependence of the branching ratio of fx — > ej on the second-generation right-handed neutrino 
Majorana mass M„ 2 in the MSSM with the right-handed neutrinos under the assumption of the MSW large 
angle solution with Vo3i = 0. We take Vd21 = —0.42 and the mu neutrino mass as 0.004eV, as suggested 
by the MSW large angle solution. The dotted line shown in the figure is the present experimental bound. 
Other input parameters are the wino mass 130GeV, the left-handed selectron 170GeV, the tau neutrino 
mass 0.07eV, and tan/3 = 3, 10, and 30. The larger tan/3 corresponds to the upper curve. 
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ji^eyin the MSSMRN with the MSW small angle solution 
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Figure 6: Dependence of the branching ratio of fx — > ej on the second-generation right-handed neutrino 
Majorana mass M„ 2 in the MSSM with the right-handed neutrinos under the assumption of the MSW 
small angle solution with Vd31 = 0. We take Vd21 = —0.04 and the mu neutrino mass as 0.0022eV, which 
are suggested by the MSW small angle solution if the mixing comes from Vd- Other input parameters 
are the same as those in Fig. (5). The dotted line shown in the figure is the present experimental bound. 
tan/3 = 3, 10, and 30, and the larger tan/3 corresponds to the upper curve. 
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Figure 7: Dependence of the branching ratio of \i — > ej on the second-generation right-handed neutrino 
Majorana mass M„ 2 in the MSSM with the right-handed neutrinos under the assumption of the 'just so' 
solution with Vd31 = 0. We take Vd21 = —0.71 and the mu neutrino mass as 1.0 x 10~ 5 cV, as suggested 
by the 'just so' solution. Other input parameters are the same as those in Fig. (5). The dotted line shown 
in the figure is the present experimental bound, tan/3 = 3, 10, and 30, and the larger tan/3 corresponds 
to the upper curve. 
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Figure 8: Dependence of the branching ratio of \i — ► e-f on the third-generation right-handed neutrino 
Majorana mass M V3 and Vd31 in the MSSM with the right-handed neutrinos. Here the tau neutrino 
mass is 0.07eV and Vd32 = —0.71, as suggested by the atmospheric neutrino result. We neglect /„ 2 here. 
The curves mean the contours on which the branching ratio of \i — > is 10 -21 , 10 -18 , 10~ 15 , and 10~ 12 , 
respectively. The shaded region is already excluded experimentally, tan/3 is set to be 3. The wino mass 
is 130GeV, the left-handed selectron mass 170GeV, and the Higgsino mass parameter \x positive. 
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Figure 9: Dependence of the branching ratio of \l — > e-f on the third-generation right-handed neutrino 
Majorana mass M V3 and Vd31 in the MSSM with the right-handed neutrinos. The input parameters 
are the same as those in Fig. (8) except that we take tan/3 = 30 here. The curves mean the contours 
on which the branching ratio of \i — > is 10~ 18 , 10 -15 , and 10~ 12 , respectively. The shaded region is 
already excluded experimentally. 
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Figure 11: Dependence of the branching ratio of r — > (jry on the third-generation right-handed neutrino 
Majorana mass M„ 3 in the SU(5) SUSY GUT with the right-handed neutrinos. Here the tau neutrino 
mass m„ T is 0.07eV and Vd32 = —0.71, as suggested by the atmospheric neutrino result. We take the 
bino mass Mi 65GeV, the right-handed selectron mass mg R 160GcV. The three curves correspond to the 
case where tan/3 = 3, 10, and 30, respectively. The branching ratio becomes larger for larger tan/3 value. 
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Figure 13: Dependence of the branching ratio of [i — ► on the third-generation right-handed neutrino 
Majorana mass M„ 3 and Vd3i in the SU(5) SUSY GUT with the right-handed neutrinos. Here we take 
the tau neutrino mass m„ T 0.07eV and the mu neutrino mass is neglected. The curves mean the contours 
on which the branching ratio of fi — > e-f is 10~ 14 , 10~ 13 , 10~ 12 , and 10 -11 , respectively. The shaded region 
is already excluded experimentally. We take the bino mass Mi 65GeV, the right-handed selectron mass 
mg R 160GeV, and tan/3 = 3. 
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Figure 14: Dependence of the branching ratio of \i — ► e-f on the third-generation right-handed neutrino 
Majorana mass M V3 and Vd31 in the SU(5) SUSY GUT with the right-handed neutrinos. The curves 
mean the contours on which the branching ratio of fi — > is 1CU 13 ,10~ 12 , and 10 -11 , respectively. 
The shaded region is already excluded experimentally. The input parameters are the same as those in 
Fig. (13) except that in this figure tan/3 = 30. 



51 




M N (GeV) 

Figure 15: Dependence of the branching ratio of /i — > e-f on the typical right-handed neutrino Majorana 
mass M N and V D 3i in the SU(5) SUSY GUT with the right-handed neutrinos. We assume the MSW 
large angle solution, which suggests m„ to be 0.004cV and Vbix = —0.42. We take the tau neutrino 
mass ra„ T 0.07eV and Vd32 = —0.60, as suggested by the atmospheric neutrino result. We assume the 
universality of the right-handed Majorana masses M v \ = M v i = M v ^ = Mjv, for simplicity. The curves 
mean the contours on which the branching ratio of fj, — ► ey is 10~ 13 , 10 -12 , and 10 -11 , respectively. The 
shaded region is already excluded experimentally. We take the bino mass Mi 65GeV, the right-handed 
selectron mass mg R 160GcV. In this figure tan/3 = 3. 
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Figure 16: Dependence of the branching ratio of \i — > on the typical right-handed neutrino Majorana 
mass M N and V D31 in the SU(5) SUSY GUT with the right-handed neutrinos. We assume the MSW 
large angle solution and the atmospheric neutrino result. All the input parameters are the same as those 
in Fig. (15) except that we take tan/3 = 30 in this figure. The curve means the contour on which the 
branching ratio of fi — ► e-f is 10~ n . The shaded region is already excluded experimentally. 
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Figure 17: Assignment of the momenta to the external leptons and the external photon in a lepton flavor 
violating diagram. An anti-charged lepton ef going into the left vertex with momentum p is annihilated 
there, and a photon with an outgoing momentum q and an anti-charged lepton e+ with an outgoing 
momentum p — q are emitted. 
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Figure 18: Patterns of the chirality flips in the lepton flavor violating diagrams in the e~l — > e^7 decay 
(i > j) . In the diagrams (a) and (b) the lepton chirality is flipped on the external lines, while in (c) and 
(d) it is flipped at a vertex of lepton-slepton-neutralino (-chargino). In (e) it flips on the internal slepton 
line. Chirality flip on the internal slepton line does not occur in the diagram with a virtual chargino 
because of the absence of the right-handed sneutrino at the low energy region. 
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